Bunching of Bell states 
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The bunching of two single photons on a beam-splitter is a fundamental quantum effect, first 
observed by Hong, Ou and Mandel. ft is a unique interference effect that relies only on the photons' 
indistinguishability and not on their relative phase. We generalize this effect by demonstrating the 
bunching of two Bell states, created in two passes of a nonlinear crystal, each composed of two 
photons. When the two Bell states are indistinguishable, phase insensitive destructive interference 
prevents the outcome of four-fold coincidence between the four spatial-polarization modes. For 
certain combinations of the two Bell states, we demonstrate the opposite effect of anti-bunching. 
We relate this result to the number of distinguishable modes in parametric down-conversion. 
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The bunching of photons is a manifestation of their 
bosonic nature. Although resulting from quantum inter- 
ference, it is insensitive to the photons' relative phase. 
Bunching of two photons on a beam-splitter was first 
demonstrated in the famous experiment by Hong, Ou and 
Mandel (HOM)Q. Later, bunched photon states were 
shown to be useful for quantum limited interferometric 
measurements and for beating the classical diffraction 
limit @. The bunching of three photons was also experi- 
mentally achieved Photon bunching became a funda- 
mental tool in quantum optics. It is used in Bell state 
analysis^, teleportation 6] and more. Recently, the ba- 
sic principles of bunching have also been found to be 
useful for quantum information processing with linear 
optics 0- An example is the nonlinear phase-shift on a 
beam-splitter for the quantum controlled-NOT gate[j|. 

Photons can bunch also in circumstances where there 
are no beam-splitters. Using the equivalence between the 
operation of beam-splitters on two spatial modes and the 
operation of waveplates on two polarization modes 0, a 
horizontally polarized and a vertically polarized photon 
can bunch when their polarization basis is rotated. Po- 
larization interference experiments are simpler as they 
involve a smaller number of spatial modes. 

The most common source for photon bunching exper- 
iments, as well as for photon entanglement, is paramet- 
ric down-conversion (PDC). As an intense pump beam 
is passing through a crystal possessing nonlinearity, 
some of its photons can split into two. For type-II non- 
collinear PDC, these two photons are emitted into two 
spatial modes (referred to as modes a and 6), and can ex- 
hibit polarization entanglement fl(i| . It is helpful to use 
pulsed pump sources due to their well defined timing. 
Recently, many experiments have used configurations 
where a pump pulse passes through the nonlinear crystal 
twice 0, EJ- In this Letter we show how four photons 
that originate from such two passes avoid being equally 
distributed between the four possible spatial-polarization 
modes by phase- insensitive destructive interference ^2] ■ 
This result has many similarities to the HOM bunching. 



Its main difference though is that instead of the bunching 
of two single photons in Fock states, the bunching occurs 
between two composite states - the Bell states. As op- 
posed to previous double-pass experiments, the result is 
unaffected by the amplitudes of two-pair emission in one 
of the two passes. We also show how, unlike in the HOM 
case of single photon states, for certain combinations of 
Bell states the bunching transforms into anti-bunching. 

Non-collinear type-II parametric down-conversion cre- 
ates the followin g bi -partite state in a single pass of the 
nonlinear crvstal|13j 
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where \m, n)i represents m horizontally and n vertically 
polarized photons in mode i. The magnitude of the inter- 
action parameter r depends on the nonlinear coefficient 
of the crystal, its length and the intensity of the pump 
pulse. The one-pair term (n=l) is the familiar Bell 
state^i). We concentrate on the case when two indistin- 
guishable photon pairs are produced (n—2). This term 
contains three elements and it is written in the above 
formalism as 

|V 2 -> = -^(|2,0)„|0, 2) 6 - |1, 1) |1, 1)6 + |0, 2>„|2,0>6) • 

(2) 

Previously, it has been shown how a four-fold coinci- 
dence between the four modes a/j, a v , bh and b v is forbid- 
den for the state of Eq. [3 when the photons of mode a 
are rotated to an orthogonal polarization basis compared 
to the photons of mode b [la flil Il7j . For example, if 
mode a is observed at the linear horizontal-vertical (hv) 
basis and mode b at the linear plus-minus (pm) 45° basis 
(or the right-left circular basis rl), an a,ha v coincidence 
collapses the state to the middle term of Eq. [21 and the 
polarization rotation bunches the b mode photons such 
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FIG. 1: The experimental setup. The pump pulses pass twice 
through the nonlinear BBO crystal, with a controlled delay 
between the passes. Down-converted photons from the first 
pass are re-injected into the crystal together with the pump 
second pass. Two BBO crystals are used to compensate for 
temporal walk-off. Photons are coupled through narrow band- 
pass filters into single-mode fibers and detected by APDs. 



that bhb v — >bp — b 2 m . Thus, a four- fold coincidence of the 
form aha v bpb m can never occur. The same result applies 
to any two orthogonal polarization bases due to the ro- 
tational symmetry of Eq. Nevertheless, when more 
than one mode are collected for each spatial-polarization 
mode (as in a double-pass configuration), these four- fold 
coincidences can be observed. 

In the experiment , we used a PDC source in a double- 
pass configuration|llj. The nonlinear crystal is pumped 
by 200 fs pulses of 390 nm wavelength from a doubled 
Ti:Sapphire laser at a repetition rate of 80 MHz. Down- 
converted photons from the first pass are redirected into 
the crystal, to meet again with the same pump pulse. The 
use of a double-pass configuration can introduce distin- 
guishability in a controlled way. As the pump back re- 
flection mirror (See Fig. ^| is translated, the timing At 
between down-converted photons from the two passes is 
changed. This timing is an extra distinguishing quantum 
number that can be varied continuously between turned 
on (by a delay longer than the coherence length t c ) and 
off (at zero delay). 

For a pure ip^ state, no four-fold events can be de- 
tected at orthogonal polarization bases and zero delay. 
Nevertheless, when At ^> t c , the temporal delay adds dis- 
tinguishability between the Bell states from the first and 
the second passes, doubles the number of collected modes 
and gives rise to four-fold coincidences. To demonstrate 
this, we define a ladder operator = (a h bl — fl^&j,)/ y/2, 
composed of creation operators of the four PDC modes. 
The operator creates a t/j^ state when applied once 
to the vacuum and a ip% state when applied twice. The 
four-photon state resulting from two passes of a pump 
pulse in the down-conversion crystal with a phase differ- 
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FIG. 2: Measured four-fold coincidence rate (squares) as a 
function of the delay At between two ip~ states. Mode a is 
detected at the hv basis and mode b is at pm. Every point is 
integrated over 25 minutes. The curve for the best fit value of 
a — 0.8 is solid while predictions for other a values are shown 
as dashed lines, inset: a fine scan around zero delay reveals 
oscillations. 

ence of u>At between them is 0] 

life) - \{4 + e iuM Ll l ) a \vac) (3) 

= - a (L\ 2 + 2e^ A * J L I t L I t I + e 2iuAt L\?)\vac) . 

The roman digits designate a distinguishing quantum 
number. The distinction can be one of many options, 
such as the photon wavelength, spatial mode or timing, 
as in this double-pass case. The terms that result from 
the first and last operators in Eq. [31 do not contribute to 
the ah<i v bpb m four-fold coincidence as they each create a 
1P2 state. Keeping the pass number labels and discard- 
ing normalization, the resulting evenly populated state 
terms, that would yield a four-fold event when observing 
modes a and b at orthogonal polarization bases are 

Infold) = (4,i4,ii 6 J>,ii & L,i + a h,n a l,i b t,i b L,n ( 4 ) 

-4,l4,II 6 p,I 6 m,II - a h,II a l,I b t,n b ln,l)\ Vac ) ■ 

Because these terms originate only from the middle term 
lJlJj of Eq. [31 their amplitudes are neither sensitive to 
the phase wAt, nor to the amplitude balance between 
the two passes. Like in the HOM bunching experiment, 
when the two input states arc indistinguishable, the pass 
indices are omitted and the amplitude for ipt-foid disap- 
pears. However, the absence of interference in the dis- 
tinguishable case revives the four-fold amplitude. As the 
delay At is scanned the distinguishability is varied, re- 
sulting in a dip for the four-fold counts, centered at zero 
delay. 

Figure [21 shows the four-fold coincidence rate for or- 
thogonal polarization bases as a function of the delay 
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between two ip^ states. As predicted, the coincidence 
rate has a dip, centered at the zero delay point. The dip 
width corresponds to the coherence time of the two ip^ 
states 0. 

Previously, the quality of the two-pairs state was de- 
fined by a, the probability for having the ibj state, as 
opposed to two distinguishable ipi states [la. Il9j : 

|V) = Va\ip2) + \/l - alVfr) ® IV-ui) ■ (5) 

This model is problematic if we assume that it originated 
from two modes as in Eq. [3J First, such a state has at 
least 60% content of ^J, i.e. a has a non-zero minimal 
value. Second, the ip^ content is the sum of two terms, 
ip2i an d V'sTn f rom the two modes (passes). Furthermore, 
two or more distinguishable modes are possible even for 
a single pass. We would like to have a model with arbi- 
trary number of modes that will support results where a 
has lower values and explain the experimental imperfect 
interference at the center of the dip. The multi-mode 
Hamiltonian and the four-photon component of the state 
produced are 

H = ^r/" /.; ■ h.C, (6) 
1 1 " d 

\Mn*)) = -^(-S^L^Iwoc), (7) 

3=1 

where k is a coupling constant that depends on the non- 
linearity of the crystal and the intensity of the pump 
pulse, and C is the proper four-photon component nor- 
malization. The distribution of Cj (real numbers) de- 
termine the number of modes involved and their relative 
weight such that their average is (c) = 1. As the operator 
U has two terms, the state of Eq. 0has 2n 2 d + rid non- 
interfering terms (3n<j originate from (i]) 2 terms and 
2nd(n,d — 1) from IAlJ k ), of which only n 2 , — nd can give 

rise to the relevant four- fold events (half of the LjL,\ 
terms). For the simple case of equally weighted modes, 
the intensity of \(j>2) is (2nd + l)/( n ci)i which is higher 
for smaller nd values due to better stimulation. After ro- 
tation to orthogonal polarization bases, four-fold events 
result only half the time and for the equally weighted 
modes case their probability per pump pulse scales as 

AMcx^i. (8) 

When only one PDC mode is collected there are no such 
events, and when the mode count increases, the rate 
peaks for nd — 2 and decays for larger numbers. When 
the time delay At is introduced in the experiment, the 
mode number is doubled such that each term from the 
first pass corresponds to a delayed term from the second 
pass. 

H = l(Hi + e™ At H n ) (9) 
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FIG. 3: Four-fold (squares) visibility and its fit (solid lines). 
Mode a is detected in the hv basis and mode b is rotated 
between hv and pm. 

At large delays At 3> t c the two corresponding modes 
are distinguishable and Eq. [5] still holds with nd replaced 
by 2nd- This case is equivalent to doubling t c (see Ref. 
|l9|). When the delay is not large enough for distin- 
guishability, the two passes of Eq. [5] interfere and create 
oscillations between zero and an upper bound defined by 
Eq. IS] Therefore, in the single-mode case, a scan will 
change the mode number between 2 and 1, thus the four- 
fold event count would present a phase-independent dip. 
When more modes are collected, the rate would oscillate 
between zero and a peak envelope. The insert in Fig. [21 
is a fine scan that reveals the expected oscillations. 

Equation [S] is a parameterization of Eq. [7|for the case 
of nd = 2 and different weights, where a — yj^(c{ + C V)- 
Thus, the four- fold rate result of Fig. [3 is composed of a 
dip contribution from a single mode per pass element and 
an oscillating contribution from a two modes per pass el- 
ement. Because the dip result of Fig. |3was integrated 
over a long time, phase fluctuations from experimental in- 
stability averaged the oscillating term. Considering this 
model with the averaging effect, the best fit for the dip 
data corresponds to a = 0.80 ± 0.05. The dashed curves 
in Fig. [5] demonstrate the predictions for different a val- 
ues. The dip is sensitive to changes in a and transforms 
into a peak for a < 0.64. 

Great care has to be taken to ensure that the col- 
lected photons originated from a single mode. Such mea- 
sures include wavelength and spatial filtering and com- 
pensation for temporal and spatial birefringence walk- 
offs. Chirped pump pulses are also a source for distin- 
guishability as their duration is longer than their coher- 
ence length 19]. Originally, it was suggested to measure a 
through the four-fold visibility, i.e. the contrast between 
four-fold coincidence count rates at similar and orthogo- 
nal polarization bases 0|. Previously published values of 
a are 37% and 83%, in Refs. 0] and 0|; respectively. 
The results of the four-fold visibility of our setup are pre- 
sented in Fig. |21 From the data we calculated a = 86%. 
The agreement with the fit value of a = 80 ± 5% from 



4 




-600 



-400 



-200 

Delay [fs] 



200 



400 



600 



FIG. 4: Four-fold coincidence rate as a function of the delay 
At between a ip~ and a (f> + states. Polarization bases and all 
other parameters are the same as in Fig. [5] Solid line is a 
Gaussian fit. Anti-bunching by a factor of 2.6 is observed. 



the dip data (Fig. |2J) and the high a value, suggests the 
validity of our two modes per pass analysis and the prac- 
tical absence of PDC states with three or more modes. 

The pure y>J" content a is also bounded from below 
for finite mode counts higher than two. The minimal 
value is a m i n = 2>/(2rid + 1). Thus, for example, a value 
of 37% indicates the necessary presence of at least four 
modes. As visibility measurements are usually affected 
by many factors, deriving a from the outcome is not suf- 
ficient for obtaining the number of PDC modes collected 
in a specific setup. The sensitive bunching measurement 
can help to remove this ambiguity. 

Bunching is not restricted to two ip^ states but can 
also occur between any two identical Bell states. Fur- 
thermore, for 8 out of the 10 possible combinations be- 
tween two of the four Bell states, bunching occurs. The 
remaining two combinations display another interesting 
result - Bell state anti-bunching. If we replace one of the 
two operators in Eq. Elby the operator K\ defined as 
the creator of a state, the last two minus signs of Eq. 
0]change into plus signs. As the delay At approaches zero 
and the four elements become indistinguishable, the four- 
fold count rate increases four times through constructive 
interference, insensitive to the phase delay as before. The 
results of such a measurement are presented in Fig. ^ 
In Table [fl we summarize the results for two polariza- 
tion cases. For the third polarization option of hvrl , 
anti-bunching occurs for ip~ ' <fi~ and for 

Bell state bunching and anti-bunching are unique 
quantum optics effects. They result from interference 
between two composite states of photon pairs. Neverthe- 
less, they are insensitive to the phase between the pairs. 
In addition, the interference is perfect even if the four 
down-conversion mode amplitudes are imbalanced. The 
Bell state bunching effects are insensitive to the ampli- 
tudes of two pairs from one pass (first and last terms of 
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TABLE I: Bunching and anti-bunching of Bell states for two 
polarization options of the modes (a, b). The letters B and A 
mark bunching and anti-bunching combinations, respectively. 



Eq. |3J), unlike any previous experiment. 

In conclusion, we presented the bunching and anti- 
bunching of two pairs of photons, each in a Bell state. 
Only the amplitudes of a single pair from each of the two 
crystal passes of a pump pulse contribute to this effect. 
It was shown how the bunching arises from varying the 
number of distinguishable modes by adding a time delay 
between the two pairs. The bunching contrast was re- 
lated to the ip2 content of the four photon state. Just 
like the HOM effect plays a crucial role in projective mea- 
surements and operations between single photons, the 
bunching of Bell states might prove to be a useful tool in 
manipulating and projecting multi-photon states. 
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